Coupling between Edge and Bulk in Strong-Field Quantum Dots 
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The maximum-density-droplet (MDD) state of quantum-dot electrons becomes unstable at strong 
magnetic fields to the addition of interior holes. Using exact diagonalization, we demonstrate that 
the first hole is located at the center of the dot when the number of electrons A'' is smaller than 
~ 14 and is located away from the center for larger dots. The separation between field strengths at 
which additional holes are introduced becomes small for large dots, explaining recent observations 
of a rapid increase in dot area when the magnetic field is increased beyond the MDD stability 
limit. We comment on correlations between interior hole and collective edge fluctuations, and on 
the implications of these correlations for edge excitation models in bulk systems. 
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The fabrication of quantum dots has now reached an 
advanced state in which the shape of the confinement 
potential and the number of electrons in a dot can be 
tuned precisely |Q. In a strong magnetic field, quan- 
tum dots display a number of interesting many body 
effects [^^, closely related to those exhibited by bulk 
two-dimensional (2D) electron systems in the quantum 
Hall regime. Correlation effects become strong because 
the single-particle level spacing in this limit is small com- 
pared to the typical Coulomb interaction energy scale. 
Some time ago we @ predicted that quantum dots in 
a strong magnetic field would form a maximum density 
droplet state with approximately uniform electron den- 
sity, close to that of a 2D system with a single filled 
Landau level {po = (27r^2)-i ^^ere £ = (Tic/eBf/'^ is 
the magnetic length), that would be stable over a wide 
field range. Indeed, features associated with the wide 
MDD state stability interval are prominent in addition 
potential vs. magnetic field traces measured by us- 
ing transport resonances to track the gate voltages at 
which electrons are added to a quantum dot. The MDD 
state is a finite-size precursor and shares many attributes 
with bulk — 1 quantum Hall ferromagnet states. 
With increasing magnetic field, po increases and Coulomb 
interactions eventually make the MDD state unstable, 
favoring a state in which the average electron density is 
lowered. Reimann et al. have concluded [ p2[ , on the basis 
of spin-density-functional calculations, that in large dots 
the density is first lowered by an edge reconstruction sim- 
ilar to those that occur in bulk systems at integer [ p^p4| 
and fractional ||l^ filling factors. In their self-consistent 
field calculations the MDD state evolves into a state with 
a modulation of charge density along the edge and a ring 
of electron density that break off from the uniform den- 
sity droplet. For smaller dots, on the other hand, Mac- 
Donald et al. have shown |^ that the average density is 
lowered by introducing holes close to the center of the 
droplet. 



This Rapid Communication is motivated by recent ex- 
periments of Oosterkamp et al., who have measured 
addition potential vs. magnetic field traces for quantum 
dots containing to 40 electrons. These authors find 
that for N larger than about 15, the area occupied by 
the quantum dot electron cloud increases rapidly on the 
high-field side of the MDD stability region. We have in- 
vestigated the quantum dot phase diagram by perform- 
ing exact diagonalization calculations, assuming Zeeman 
fields large enough to ensure full spin polarization. We 
find that holes are added to the interior of the dots with 
increasing field, increasing the electron cloud area; for 
N > Nc {Nc « 13 — 14) the interior holes stay away 
from the center of the dot, while for N < Nc the holes 
form a puddle at the center of the dot. The phase dia- 
gram we obtain for hole number vs. field, is in qualitative 
agreement with experimental data. We find that the rate 
at which holes are added with field increases for larger 
dots, explaining the observations of Oosterkamp et al.. 
We also find that correlations develop between the mo- 
tion of holes in the interior of the dot and collective edge 
excitations of the dot, and discuss the nature of the corre- 
sponding coupling in bulk integer filling factor quantum 
Hall systems. 

We consider a system of 2D electrons confined by a 
parabolic potential, V{r) — m*il^r^ /2. and confine 
our attention here to the strong magnetic field limit, 
fl/iVc <S^ 1. (wc = eB/m*c, where B is the magnetic 
field perpendicular to the 2D layer.) In this limit the 
symmetric gauge single-particle eigenstates in the lowest 
Landau level [ p^ are conveniently classified by an angu- 
lar momentum index to > 0: 

e,„ = ?ici;c/2 + 7(m -I- 1), (1) 

where 7 = mil^P We will drop the constant kinetic en- 
ergy from subsequent discussion. The Hamiltonian is in- 
variant under spatial rotations about the axis that is per- 
pendicular to the 2D plane and passes through the center 



of the dot. It follows that the total angular momen- 
tum Mz is a good quantum number. The Hamiltonian of 
the dot in the fully spin-polarized, lowest-Landau-level, 
Hilbert space is 



77: 
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The characteristic energy for the two-particle Coulomb 
matrix elements is e^/e£, where e is the host semiconduc- 
tor dielectric constant. The eigenstates of this Hamilto- 
nian depend only on the dimensionless ratio between the 
two competing terms, 7 = ^/{e^/et). The MDD state 
is the eigenstate of this Hamiltonian with the minimum 
angular momentum allowed by the Pauli exclusion prin- 
ciple, Mz -> M^^^ = N{N - l)/2, and is the ground 
state of this Hamiltonian for strong confinement (large 
7) . We find it useful below to denote occupation number 
eigenstates by lists ordered by increasing angular mo- 
mentum. For example, in this notation the TV = 4 MDD 
state is (•, •, •, •, o, o, ...), where the symbols • and o de- 
note occupied and unoccupied single-particle states. 




plots are, for each N, near the MDD stability limit. In 
HF theory, the MDD becomes unstable when an unoccu- 
pied interior quasiparticle energy exceeds the edge quasi- 
particle energy, i.e. when eHp{m*) > enpiN — 1) for 
m* < N — 1. We see from Fig. 1 that holes are first 
introduced near the center of the dot < 12 and gradu- 
ally move outward as N increases. This qualitative pre- 
diction of Hartree-Fock theory is generally confirmed by 
the exact diagonalization calculations we discuss below. 

To understand the role played by correlations we first 
discuss a small dot for which an accurate analytic calcu- 
lation is possible. For TV = 3 the MDD has M*^^^ = 3 
and HF theory predicts that when 7 is reduced the first 
hole is introduced in the m* = state, increasing the 
iV = 3 ground state angular momentum to Mz = 6. The 
exact single hole ground state is obtained by diagonaliz- 
ing the many-particle Hamiltonian in the space of states 
with Mz = 6. Fig. 1 suggests that the probability of oc- 
cupying single-particle states with large m is small, since 
quasiparticle energies increase rapidly outside the MDD. 
If we can neglect the possibility of occupying states with 
m > nic — 4, the only additional state in the Hilbert 
space with Mz = 6 has the interior hole at m = 1 rather 
than TO = and excites an edge magnetoplasmon at 
the edge |^,|l8| by making a particle-hole excitation from 
TO = 3 to TO = 4. The many-particle Hamiltonian matrix 
in this two-dimensional Hilbert space is 
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where 
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Umi7n2ni3m4 = < TOiTO2|y|TO3TO4 > 

— < TOiTO2|y|TO4TO3 > . (4) 



FIG. 1. Hartree-Fock quasi-particle energies eHFim) 
of the MDD state as a function of m for different 
values of A'^. 7 is chosen, for each N, so that the 
MDD is close to its instability point; 7 =0.077, 0.08, 
0.082, 0.085, 0.0875, 0.0925, 0.0975, 0.1, 0.108, 0.113 
for = 30, 28, 26, 14, 12, respectively. The inset dis- 
plays the HF phase boundary between the MDD state 
and one- hole states for a dot with HQ = 3meV. The di- 
mensionless confinement strength parameter 7, the field 
in Tesla B, and the confinement oscillator frequency in 
meV are related by B[Tesla] = [0.131(hQ[meY]f /jf^^ . 

In Fig. 1 we have plotted Hartree-Fock quasiparti- 
cle energies ,eHFim), for to = to to = 40 for MDD 
states with N = 12, 14, 28, 30. {eHF{m) = 7(to + 

1) + J2Z~=o Um,m' where Um,m' =< mm'\V\mm' > - < 
mm'\V\m'm >) The confinement strengths used in these 



The interaction contributions to the matrix elements in 
Eq. 3 in e^/ei units are 0.940 and 1.114 respectively 
along the diagonal and 0.083 in the off-diagonal term. 
The off diagonal term introduces correlation between hole 
motion and collective edge excitations in this small dot, 
reducing the ground state energy to 97 -|- 0.907e^/e^. 
The iV = 3 ground state energy, in this approxima- 
tion, is 97 + 0.907e^/e£, compared to the MDD state 
energy which is 67 -I- 1.204e^/e£. Note that correlations 
between bulk and edge shift the MDD stability limit 
from 7 = 0.088 in the Hartree-Fock approximation to 
7 = 0.099. An exact diagonalization for TV = 3 confirms 
the accuracy of this Hilbert space truncation. 
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FIG. 2. Occupation number distributions for the 
one-hole MDD states of A*' = 6, 12, 13, 14, 20, 30 quantum 
dots. 



FIG. 3. Occupation number distributions for one-hole 
(circle), two- hole (square), and three-hole (diamond) 
MDD states for iV = 9 and iV = 17. 



Fig. 2 summarizes our numerical exact diagonaliza- 
tion results for one hole states at larger N by plotting 
mean-occupation numbers for all single-particle states. 
For iV = 6, the HF approximation places the hole at 
TO* — 0, yielding a total = 21. The exact single- 
hole state occurs at the same value of M^, but the hole 
position and the dot edge have correlated quantum fluc- 
tuations, evidenced in Fig. 2 by the finite probability for 
finding the hole in the m = 1 state. Equal time correla- 
tion functions, not illustrated here, demonstrate that the 
hole has a higher probability of being on the same side of 
the system as edge density modulations. For N = 12, we 
again find that the exact single-hole state has the same 
Mz as the Hartree-Fock state with m* — 0, — 78, 
and correlated bulk-hole edge excitations. This pattern 
continues for larger particle numbers with the nominal 
hole angular momenta moving out at larger N. For 
N = 30, for example, the nominal angular momentum 
of the single-hole is m* = 12. Note that the hole angu- 
lar momentum has fiuctuations only to the high-angular 
momentum side of its nominal value. This property is a 
combined consequence of angular momentum conserva- 
tion and the chiral |19| nature of the collective edge exci- 
tations with which it is correlated. The one-hole nominal 
locations for N = 6, 12, 13, 14, 20, 30 m* = 0, 0, 0, 1, 5, 12 
respectively. The total hole occupation number deter- 
mined by summing the depleted occupation numbers in- 
side the quantum dots in Fig. 2 is accurately quantized 
at 1, demonstrating that holes and edge excitations are 
correlated but that charge fluctuations from bulk to edge 
or vice versa is negligible. 



For each particle number N, the number of holes in the 
ground state increases as 7 decreases, with the number 
of bulk holes accurately conserved to larger hole numbers 
for larger N. In Fig. 3, we plot occupation number dis- 
tributions for the one-hole two-hole and three-hole states 
at iV = 9 and iV = 17. We flnd that the average number 
of interior holes in the nominal two and three hole states 
deviates from integer values by 0.007 and 0.05 for = 9 
and 0.003 and 0.006 for N = 17, demonstrating accu- 
rate interior hole-number conservation. The one-hole and 
two-hole Mz values are 45 and 50 for iV = 9 and 148 and 
158 for N = 17, giving nominal flrst-hole and second-hole 
TO* values of and 4 for = 9 and 5 and 7 for N — 17. 
It is clear from Fig. 3, however, that correlations within 
the two-hole system make these assignments, based on 
the sequence of ground state angular momenta values, 
less meaningful. For States with two or more holes, cor- 
relations among the holes and correlations between the 
holes and the edge system are both important. 
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FIG. 4. Phase diagram near the stabihty hmit of the 
MDD state. The symbol MDD stands for the maximum 
density droplet, while IHC and IHR represent represent 
respectively centered (m* — 0) and ring (m* > 0) sin- 
gle-hole states. 2H and 3H represent states with two and 
three holes. States with more than three holes have not 
been studied but will occur at larger fields. Classifying 
states near the MDD limit by the number of holes is mean- 
ingful when charge fluctuations between the interior of the 
MDD and its edge are negligible. 

To facilitate comparison with experiment, we have 
plotted a phase diagram of MDD related states in Fig. 4 
vs. particle number N and field, parameters that are 
known experimentally, using a typical value for the 
parabolic confinement strength hfl — "imeV . We have 
not considered states with more than three holes because 
of the growing difficulty of carrying out accurate exact di- 
agonalization calculations, and expect that the region at 
the top right of this phase diagram should be occupied by 
states with four and more interior holes. Note that the 
spacing in field between hole number addition points be- 
comes small for large dots, explaining the sudden increase 
in quantum dot area observed by Oosterkamp et at. The 
stability regime of multiple hole states is seriously under- 
estimated by Hartree-Fock and density-functional-theory 
calculations which are unable to account for the strong 
correlations possible in the quantum Hall regime. 

The sign of cusps in the magnetic field dependence 
of the addition spectrum contains valuable information 
about the shape of the phase boundaries. Oosterkamp et 
al measured is the magnetic field dependence of the 
addition spectrum given by /ijv = — En^i + huJc/2, 
where En = 7(iV + M^) + U{N,M^) is the groundstate 
energy of an TV-electron system. (Here U{N,Mz) is the 
total interaction energy). As the applied magnetic field 
increases, groundstate level crossings in the — 1 and 
N particle systems lead to positive and negative cusps 

, respectively. For > 8 the phase boundary between 
the MDD and one-hole states is an increasing function of 
magnetic field (see Fig. 4). This implies that the A^ — 1- 
electron dot will cross the phase boundary before the 
A'^-electron dot does. The measured upward cusps are 
consistent with this expectation. Upward cusps are also 
seen between one-hole and two-hole states. 

The strong correlations that we find between interior 
holes and the edge excitations of quantum dots also 
applies to bulk quantum Hall systems. These corre- 
lations are important because of the long-range of the 
Coulomb interaction, which causes charge fluctuations in 
the two subsystems to interact strongly. Recent scan- 
ning probe studies of 2D electron systems in the quan- 
tum Hall regime have made it clear that real experi- 
mental samples always have low-energy interior hole and 
electron low-energy degrees of freedom. At low temper- 
atures charge carriers in the bulk are localized so these 



disorder induced degrees of freedom do not influence the 
quantum Hall effect. They may, however, influence other 
properties of the chiral edge state system. Collective ex- 
citations at the edge, will induce electric potentials in 
the interior that will excite the localized hole system. 
Similar effects cause semiclassical edge magnetoplasmon 
pl[ excitations to be influenced by the bulk conductivity. 
The possible importance, at experimental temperature 
and voltage scales, of these interactions for the tunneling 
density-of-states Q and other properties measured at 
the edge deserves more careful consideration. 
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